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Abstract 

The azimuthal anisotropy of high pr particle production in heavy ion collisions is a sensitive probe of the jet 
quenching mechanism. Recent PHENIX measurements for Au+Au collisions at ^JsNN = 200 GeV show a V2 
signal that exceeds the pQCD energy loss calculations up to ~ 10 GeV/c, challenging the traditional perturbative 
picture of the energy loss process. Here, we present an update and details of that measurement, as well as new high 
Pt measurements at ^/snn = 62 and 39 GeV. These measurements not only provide an important constraint for 
understanding the path-length dependence of jet energy loss and the role of initial collision geometry, but also allow 
a search for the onset of jet quenching as -sfs^N is varied. 



The discovery of the suppression of both high pr single particle yields [1] and back-to-back jet yields in two 
particle correlation measurements [2] ("jet-quenching") is considered one of the most convincing piece of evidence 
that de-confined matter (or a quark gluon plasma [QGP]) is created in relativistic heavy ion collisions at RHIC. The 
study of jet quenching has developed into a full-fledged field with precise diff'erential measurements and quantitative 
theoretical calculations. However, our current understanding of jet quenching faces several challenges: on the one 
hand, the AdS/CFT approach to describe jet-medium interactions seems to be favored by several measurements (e.g. 
Raa for non-photonic electrons [3] and high px V2 [4]) over the traditional pQCD based treatment; on the other hand, 
diff"erent models which operate inside the pQCD framework, show large diff'erences in their predicted quenching 
parameters, even though they are embedded into the same evolving medium characterized by 3D hydro and tuned to 
reproduce central Raa [5]. Precisely measured anisotropy at high pr, can shed new insights on both challenges. 

In 2007, PHENIX accumulated a large dataset of Au+Au Collisions at ^fs^N = 200GeV with several Reaction 
Plane (RP) detectors which spanned a broad pseudorapidity (77) range. They included an inner ring (RXNin) (1.0 < 
77 < 1.5) and an outer ring (RXNout) (1.5 < 77 < 2.8) of a newly installed reaction plane detector (RXN), a beam 
beam counter (BBC) and the muon piston calorimeter (MPC) (3.1 < 7/ < 3.9). These RP detectors not only provided 
precise measurements of the Event Plane (EP) but also enabled a systematic examination of non-flow eff'ects, such 
as a possible jet bias. In 2010, RHIC initiated the energy scan program and PHENIX collected several datasets for 
Au+Au collisions at ^JsNN = 62 and 39 GeV with significant statistics with the same set of RP detectors as in 2007. 
The Measurements of V2 from all of these datasets are discussed in the following. 

The techniques used for data analysis are adopted from Ref. [4]. Briefly, photons are identified and measured 
using the electromagnetic calorimeter (EMCal). The invariant mass for photon pairs is sampled over events and a 
statistical subtraction is implemented to remove the combinatoric background; the shape of the latter is estimated 
via the widely used event-mixing technique. An integral of the subtracted invariant mass distribution gives the 
yield. The emission angle of relative to the EP (A0) is folded into [0, 7t/2] (due to the symmetry in the event- 
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averaged collision geometry; for in-plane and nil for out-plane) and further divided into six angular bins; yield is 
measured individually in each of the angular bins. The resulting modulation of ''jfi yield vs. A0" is Fourier analyzed 
[6] to obtain the V2 term. A fit gives the so-called v^"^ and the true corrected V2 is the quotient v^^^ Ictrp, where cr^^p's 
are RP resolution correction factors. These factors characterize the dispersion of the measured EP (^) relative to the 
"true" RP and are defined as ctrp = <cos 2Q¥ - ^^rp)) [6]. 

The t:^ V2 results for the three collision energies are shown in the left panel of Fig-1. For both central (0 - 20%) 
and mid-central (20 - 40%) events the V2 values for the respective beam energies are quite similar. However, only the 
200 GeV data has the statistical significance to probe higher pr's (> 6.0GeV) where jet fragments dominate hadron 
production. Below 6 GeV V2 is described by hydrodynamics at low pr (< 2.0GeV) and recombination at intermediate 
Pt (2.0 - 6.0GeV); the good agreement of V2's across collision energies indicates that V2 is saturated in these two 
regions. To facilitate the investigation of V2 saturation, the V2 vs. ^fsNN plot is presented on the right panel of Fig-1. 
At low Pt, our measurements, while showing good agreement with charge hadron results for 200 and 62GeV from 
[7], add one more point at 39GeV. Combined with V2 at even lower ^JsNN at 17 and 3GeV (measured by CERES and 
E895 respectively), it seems that a transition to saturation is observed. At intermediate pr, V2's are saturated down to 
39GeV as well. 
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Figure 1: Left panel: V2 vs. pr for 200, 62 and 39GeV for centrality - 20% and 20 - 40%. Right panel: a V2 energy scan with and charge 
hadron (from [7]; all charge hadron data are of pr 1.5 - 2.0GeV). Statistical errors are indicated by bars and systematic errors by square brackets. 

In the following we focus on the high pp region and discuss results from 200 GeV only. A complementary 
measurement is V2 and Raa for the r] meson, which is reconstructed from two photons as well and analyzed in the 
same manner as for the n^. The r] has a larger mass so its yield can only be extracted for pp > 4.0GeV in PHENIX; 
this is the region where jet-medium interactions play a central role, r] V2 and Raa as a function of pp are shown in 
Fig-2 for central (0 - 20%) and mid-central (20 - 60%) collisions; the results are overlaid to aid a comparison. The 
measurements show solid agreement within one standard deviation. This observation indicates that mass and quark 
contents have little effects on high pp V2 and Raa, and probably the underlying jet-quenching mechanism as well. 

Fig-3 shows t:^ V2 and Raa vs. Npan (number of participants) for 6 < pp < 9 GeV/c. The results from four popular 
pQCD model calculations are also plotted for comparison. Note that the parameters for these models are tuned to 
match Raa for the most central collisions. It is clear that while all the models successfully reproduce Raa for a good 
range of Npart \ they substantially under-predict V2. Moreover, since V2 and Raa are anti-correlated as indicated by 



Hhe WHDG model undershoots the data points because it is tuned in the pr range 5 ~ 18 GeV/c and tends to undershoot the low pr end. 
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Figure 2: 77 (black square) Raa (left) and V2 (right) vs. pr, overlaid on results (grep circle) for comparisons. Two centrality ranges are presented: 
- 20% (top) and 20 - 60% (bottom). 



the solid edge on the WHDG band on the figure, the measurements of V2 and Raa can not be simultaneously described 
by these models. This discrepancy might be rooted in a combination of diff'erent eff'ects. In Ref. [8] diff'erent initial 
collision geometries were considered; instead of the default usage of Glauber Model geometries, the CGC type of 
geometry as well as its event-by-event fluctuations were taken into account. However, the resulting enhancement 
of V2 was not able to account for the whole gap. Subsequently, a model calculation based on AdS/CFT like jet- 
medium interactions was tested. The idea is that all the models above are based on pQCD; the underestimation of V2 
may suggest that a stronger jet-medium interaction and thus, stronger path-length (L) dependence of jet-quenching 
is necessary. The AdS/CFT type of path-length dependence is L^, stronger than the pQCD counterpart L^. This 
modification is remarkably eff'ective; combined with CGC geometry and event-by-event fluctuation, the reproduced 
V2 matches data points well. 




Figure 3: V2 and Raa vs. Npart- Four pQCD models are plotted for comparison: ASW, AMY and HT from [5]; WHDG from [9]. The upper 
bound of dNg/dy in WHDG (1600) is indicated by the solid edge on the band. 
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To further investigate path-length dependence of jet-quenching, angular dependent Raa is introduced and denoted 
by Raa(^(P)- It combines the absolute suppression level Raa with path-length variations from in-plane to out-plane. 
Fig-4 shows high /?7 (7 ~ 8 GeV/c) Raa(^(P) measured for 6 A0 bins in 6 centrality bins from the most central up 
to 60%. For each A0 and centrality bin the path-length integral (m = 1,2) is calculated based on an initial geometry 
represented by p. On the right panel of the figure, where assumptions of an AdS/CFT like jet-medium interaction 
(I2) and CGC geometry with event-by-event fluctuation (p^q^) are made, Raa scales remarkably with path-length 
integrals. This scaling suggests that the inclusive Raa and V2 for diff'erent centralities (or Npart) will be automatically 
and simultaneously described. 
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Figure 4: Raa{^<P, centrality) vs. 7^: Raa{^<P, centrality) is measured in 6 A0 bins and 6 centrality ranges; Im is the path-length intergral and 
m = 1,2 corresponds to pQCD and AdS/CFT like jet-medium interactions respectively, p represents initial condition: p^^^^ for default Glauber 
geometry and p^^^ for CGC geometry with event-by-event fluctuation. 

In summary, tt^ V2's for 200, 62 and 39 GeV are measured. In the low and intermediate pr range, V2 is saturated 
across energies. For higher pr values we found that current pQCD models are challenged to reproduce the V2 and 
Raa simultaneously. However, an AdS/CFT-based calculation that encodes an path-length dependence with an 
event-by-event fluctuating CGC geometry is favored by the measurements. 
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